Class of exactly soluble models of one-dimensional spinless fermions and its 
application to the Tomonaga-Luttinger Hamiltonian with nonlinear dispersion 
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I. INTRODUCTION 

The Tomonaga-Luttinger model is the most important 
model of one-dimensional interacting electron physics. 
The model describes two Fermi points of right-moving 
(R) and left-moving (L) electrons. The electrons interact 
through a two-body short-range potential: 



Ht~L — -ffkin + -ffint , 

#kin = ivF I dx[: V>iV?/> L : -: ^r,W r : ) , 

-l-L/2 V 
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It is shown that for some special values of Hamiltonian parameters the Tomonaga-Luttinger 
model with nonlinear dispersion is unitary equivalent to the system of noninteracting fermions. For 
such parameter values the density-density propagator of the Tomonaga-Luttinger Hamiltonian with 
nonlinear dispersion can be found exactly. In a generic situation the exact solution can be used as a 
reference point around which a perturbative expansion in orders of certain irrelevant operators may 
be constructed. 



Thus, some researchers began developing 'fermionic' 
approaches to the model. In a variety of contexts the 
one-dimensional electron systems were described in terms 
of the composite fermionic degrees of freedom [3, 4, 5]. 
More to the topic of the present paper are Ref. [6, 7, 8] 
which we briefly discuss below. 

In Ref. [6, 7] the authors studied the effect of non- 
linear dispersion v' F on the density-density propagator of 
the Tomonaga-Luttinger Hamiltonian. They avoided us- 
ing the bosonization and instead worked with the bare 
fermionic degrees of freedom. 

The quasi-one-dimensional case was investigated in 
Ref. [8] where the hybrid boson-fermion representa- 
tion for the quasi-one-dimcnsional Hamiltonian was pro- 
posed. The latter representation describes the high- 
energy one-dimensional degrees of freedom in terms of the 
Tomonaga-Luttinger bosons while the low-energy three- 
dimensional degrees of freedom as fermionic quasiparti- 
cles. 

Despite their success, both approaches are not sys- 
tematic. Rather, they were custom-made for particular 
tasks at hand. The systematic method to handle the 
generalized problem (l)-(5) was proposed by the author 
in Ref. [9]. The main idea of the method is to apply 
a certain canonical transformation to (l)-(5) which kills 
marginal (in the renormalization group sense) interaction 
-Hint. Such transformation maps the original model on 
the model of fermionic quasiparticles weakly interacting 
through an irrelevant operator. Due to the irrelevance of 
the quasiparticle interaction the usual perturbation the- 
ory is applicable and the quasiparticles near the Fermi 
points are free. 

This approach allows one to calculate the retarded 
density-density correlation function of the Tomonaga- 
Luttinger model with nonlinear dispersion to the zeroth 
order in the quasiparticle interaction [9]. 

The developments described above solve the long- 
standing problem of the Tomonaga-Luttinger model with 
the nonlinear dispersion. The purpose of this paper is re- 
finements of the method: we will show that if, in addition 
to the nonlinear dispersion, one supplies the Tomonaga- 
Luttinger Hamiltonian with a certain type of irrelevant 
interactions the Hamiltonian becomes exactly soluble. 
The system at and near this solubility point is studied 



H mt = I dxdx'g(x - x')pl(x)pr(x'), 



(1) 
(2) 

(3) 
(4) 



where function g{x) vanishes for A|x| ^> 1. The param- 
eter A plays a role of the ultraviolet cutoff. The symbol 
: . . . : stands for the normal ordering with respect to the 
noninteracting ground state. 

The kinetic energy (2) of the model (1) is a linear func- 
tion of the electron momentum. It is a necessary condi- 
tion for the exact solubility of H^l- when the kinetic 
energy is linear one can apply bosonization to find the 
excitation spectrum and the correlation functions of the 
model [1, 2]. 

However, there are situations in which we must extend 
the model by including the nonlinear dispersion term: 



dx 



V^V^ R : 



(5) 



Such generalization immediately introduces serious com- 
plications: the bosons become interacting. 

Another extension of the system (1) where the 
bosonization fails is the quasi-one-dimensional conduc- 
tor: an array of one-dimensional Tomonaga-Luttinger 
wires coupled by transverse hopping and interaction. 
Bosonization works poorly for the transverse hopping has 
no simple expression in terms of the Tomonaga-Luttinger 
bosons. 
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in this paper. 

We proceed by noting that the classical formulation of 
the Tomonaga-Luttinger model with nonlinear dispersion 
(l)-(5) is somewhat inconsistent. It is trivial to check 
that the scaling dimension of H n \ is equal to 3. Indeed, 
the field operator has dimension of 1/2, the dimension 
of V is unity. Therefore, the dimension of the whole 
operator is 3. In addition to H n \, there is one and only 
one operator whose scaling dimension is also 3: 

H! nt = ig'j dx [pr, (:4(VVO: " = (W-i> L :) (6) 
-Pl(:4(Wr):-:(V4)V<r:)" • 

From the renormalization group view point if we are in- 
terested in the effect of H n \ we must include H[ nt as well. 

Initially, it looks as if such modification leads to sub- 
stantial additional difficulties. Yet, despite this expecta- 
tions, we will see that the operator H[ nt does not make 
our problem more complicated at all. Moreover, if g' , v' F , 
g and vf satisfy some special relation the Hamiltonian 

H cx = H T L + H nl + H! nt (7) 
+ (irrelevant operators with scaling dim > 4) 

is equivalent to the system of noninteracting fermions. 
Thus, the system becomes exactly soluble. 

The exact solubility gives us a chance to access qual- 
ities otherwise hidden from the analytical investigation. 
For example, the spectrum of the Hamiltonian and the 
statistics of excitations can be determined. Additionally, 
the density-density propagator can be calculated as well. 
It is proportional to the propagator of the free fermions. 

Of course, a generic model (which we denote as Hq, 
where 'G' stands for 'generic') deviates from the exactly 
soluble case. However, it can be shown that if such devi- 
ations are weak we can study the generic model with the 
help of perturbation theory in orders of (Hq — H cx ). In 
particular, the propagator of the exactly soluble Hamilto- 
nian is the zeroth order approximation to the propagator 
o{H G . 

The paper is organized as follows. The exactly solu- 
ble model is derived in Sect. II. Since the latter Section 
is very technical the reader might get lost in the details. 
For those who are not specifically interested in all the 
minutia there is Sect. Ill where we give concise nontech- 
nical overview of the presented derivations. This Section 
is fairly self-contained so that the reader ready to set- 
tle for just a heuristic argumentation may skip Sect. II. 
The situation of weak deviation from the point of ex- 
act solubility is discussed in Sect. IV. Sect.V is dedicated 
to the application of Sect. IV ideas to calculations of ob- 
servables. In Sect. VI we compare our approach with the 
results available in the literature. Technically involved 
calculations are relegated to Appendix. 



II. EXACTLY SOLUBLE MODEL 

We start with the quadratic Hamiltonian: 

Ho = -ffkin + H n i + . . 

J2 J dx{ l pvr.^(^ P y+v F -.(y^ P )(y^ P y} + 



(8) 



where ellipsis stand for quadratic terms which contain 
more space derivatives. We must assume their presence 
to avoid creating spurious Fermi points besides those lo- 
cated at q = 0. From the practical point these terms 
cause us no trouble. They will be discussed in Sect. IV, 
together with other terms whose scaling dimension is 4 
and higher. The summation in eq.(8) runs over the chi- 
rality index p which is equal 1 for left-movers and -1 for 
right-movers. 

As we have mentioned above the colons denote the nor- 
mal ordering with respect to the noninteracting ground 
state |0). Specifically, the normal order is defined by the 
equations: 

^l{x)%{y) = :ipl(x)ip p (y):+s p (x-y), (9) 
^ p (x)%(y)^ p (x')%(y') = ^l(x)^ p {y)^l(x')%{y')ilQ) 
+s p (x - y):i>l(x')i> p (y'): +s p (x' - y'):ip f p (x)ip p (y) 
+s p (x - y')-4 p {y)i>l(x'): +s p (y - x'):i>l(x)i> p (y') 
+s p (x - y)s p (x' - y') + s p (x - y')s p (y - x'), 

^(^ = (01^(^(0)10) =(11) 

Obviously, the Hamiltonian (8) is trivially soluble. The 
spectrum of the fermionic excitations is 



£ pq = pi>Fq + v F q 2 + 0{q 3 ). 



(12) 



The contribution 0(q 3 ) comes from the terms of cq. (8) 
denoted by the ellipsis. These terms guarantee that e pq — 
for q = only. 

The density-density propagator for the model (8) 



v o 



to 2 + (v F q - v' F q 2 ) 2 
4TTv' F q ° \ iv 2 + (vpq + v' F q 2 ) 2 



■log 



can be calculated without much effort. 

Now we define the unitary transformation U: 



n 



p 



u 



P Pp-qP-pqi 



Lq 

2tt' Ppq 



dxe lqx p p (x), 



and transform with its help the operator H : 
H cx = UH Ul 



(13) 

(14) 
(15) 

(16) 
(17) 
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The subscript 'ex' stands for 'exact' since the Hamilto- 
nian H CK is unitary equivalent to the quadratic Hamilto- 
nian (8) which implies exact solubility of H ex . 

The coefficients a's from eq.(15) uniquely specify the 
operator U. We impose on them the following conditions: 



-qi 



- a l\q=0 fOT M < A ' 

a q — > for \q\ > A. 



(18) 
(19) 
(20) 



These conditions mean that U when acting in the Fock 
space of Hamiltonian (8) cannot create highly excited 
single-fermion states. Loosely speaking, the action of U 
is confined to the region \q\ < A near the Fermi points. 
This guarantees that the interactions of H ex have proper 
ultraviolet cutoff of A. 

In order to find H cx explicitly we need the following 
commutation rules: 



[Ppqj Pp'~ <?'] — Spp'Sqqipriq, 



Sppf 6 



(21) 
(22) 



which are proven in [1, 10]. 

Next, let us study how operator U acts on different 
observables. With the help of cq.(21) it is possible to 
show that for q ^ 0: 



U p pq U ] = u(a q )p pq + v(a q )p 



u(a q ) = u q = ch a q , v(a q ) = 

For zero modes we have: 

UNpU^ = N p , 

N P = Ppq\ q = - 



■pqi 

sh a q , 
Z = 1. 



(23) 
(24) 
(25) 

(26) 
(27) 



Such transformation rule is a consequence of definition 
(15), in which commutes with N p . 
It is convenient to define the following symbols: 



Pp( x ) = J dyu{y)p p {x -y) = (u* p p )(x), (28) 

Pp(x) = J dyv(y)p p (x - y) = (v * p p )(x), (29) 

u(x) = L- 1 Uqe iqx - f ^ V x . ( 3 °) 
v(x)=L- iy £v q ^= J Y^ iqx - 



(31) 



Here the asterisks stand for convolutions of two functions. 

Due to eq.(20) we have u q = 1 and v q — for \q\ 3> A. 
Therefore, v(x) is well defined, whereas u{x) is singular: 

u(x) = 6(x) +w(x), (32) 

dq - ~ iax (33) 

(34) 



2tt 



-w„e 



q^ 1 



where w(x) is well-defined. 

Both v(x) and w(x) are even functions of x. They vary 
slowly on the scale of (1/A) and vanish for \x\ 3> (1/A). 

We express the action of U on p p (x) in the following 
manner: 

U Pp (x)U^ = p u p {x) + p%(x) + Spp (35) 
= P P (x) + P p (x) + P-p(x) + 5p p0 , 

Pp "=w* p p , (36) 



5 p. 



'p0 



- [(1 - u )N p - v N_p] (37) 



-— [w a N p 



v N. 



Uq = Ug|,=0, V = V q \ q=0 , W = W q \ q=0 . 



(38) 



The obtained result for UpU^ can be used to find 
UH^U^ . We prove in Appendix that: 



| [:iil(x)(V%(x)):-:Wl(x)W p (x):] (39) 
= 7T lim {p p (x)p p (y) -bp(x-y)}, 

bp(x) = {Q\pp(x)pp(0)\0)= S 2 p (x). (40) 

These formulas express the kinetic energy density in 
terms of the density operators. Since we know how the 
latter are transformed under the action of U we deter- 
mine UH K i n lI^ easily, as the derivation below demon- 
strates. 

First, introduce the notation: 

: Pp(x)p P (y)- = pp(x)p p (y) - b p (x - y), (41) 
: Pp(x)Pp(y)p P {z): = pp{x)pp(y)pp(z) - (42) 

b P ( x - v)pp( z ) - b p( x - z )pp(v) - K(y - z )pp( x )- 

In the bosonization framework the above equations define 
the normal ordering of the bosonic operators. In this pa- 
per we do not introduce the Tomonaga-Luttinger bosons. 
Thus, the symbols : pp: and : ppp: must be understood 
as a shorthand notation for expressions like eq.(41) and 
eq.(42) without further hidden field-theoretical meaning. 
The reason for using normal ordered products is that 
such objects are well defined analytically at any values 
of x, y and z while unordered products have singulari- 
ties when arguments coincide (e.g, the unordered product 
p P (x) p p (y) has a second order pole at x — y). 

The new notation allows us to reformulate eq.(39) in 
a more compact way: 



f [:$(W P ):-:(V$W P :] = 7r:# . 



(43) 



w q = u q - 1 — ► for \q\ > A, 



In the bosonization context this equality corresponds to 
the bosonization of the kinetic energy [1, 10]. 
Using eq.(39) and eq.(43) we find: 

UH kin U^ =7:i, F J2 [ dxU: p 2 p . U f = (44) 
p J 

ttv f / dx\^ lim U {p P (x)pJy) -b p (x-y)}U^. 

J V^x 
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The expression under the limit sign is equal to: 

U{p p (x)p p (y)-b p (x-y)}U* = (45) 
(pp(x) + P- P (x) + 6p p0 ) ( P ;(y) + P %(y) + S Pp0 ) 

-b p (x - y). 

In this equation we need to normal order p operators 
according to definition eq.(41): 

p P l (x)p p >(y) = :p p i (x)pZ(y):+(u*b p *u)(x-y). (46) 

It is easy to show with the help of eq.(25) that for any 
function f(x) 

(u* f * u)(x) = f(x) + (£*/* v)(x). (47) 

Therefore, eq.(46) may be written differently: 

p u p (x)p u p (y) = :p;(x)p;(y):+ (48) 
b p (x -y) + (v*b p * v)(x - y). 

Finally: 

p\{x)p\{y) = : p%{x)p%{y): +(v * b^ p *v)(x- y\A9) 
Collecting all terms together one finds: 

UH kin U<< = tt5 f f dxJ2- (Pp + P-p + <W>) 2 : (50) 
J v 

+2ncv F A 2 L, 

where c is the dimcnsionlcss non-universal c-number: 



The term h[^ t is a sum of irrelevant operators whose scal- 
ing dimensions are 4 and higher. That is why we put such 
superscript at h: 

CO ,. 

h$ = 2 E < A_2 " J dx ■ Pp(V 2n P P )-- , (55) 



A 2 " ( 

= J2nf\ J dxx (v*v)(x), (56) 
c o= f dx{v*v){x) = v 2 . (57) 



We establish: 



J dx{:{plf:+:{p^?:}= J dx(l + 2c* ): p 2 p . +h$(58) 



= (l + 2cg)^ / dx:^ p :+h£l 



Here we used eq.(43) to transform : p 2 :. 

The : plPr: term of eq.(52) can be written as: 

2nv F J dx{plpl + plpl) = (59) 

J dxdx'g(x - x / )p L (x)p R (x'), 

g(x) = 4itv-f(v * u)(x) g q — Aitv F v q u q . (60) 

Ultimately, we obtain the following expression for 
UH kiD W: 



/dq r L / 2 
Jr^j^^K- (51) UH kin U^ = iv F J dx^lViP L :-.^ R V^ R :) (61) 



In cq.(50) contributions proportional to Sp pn or (Sp p o) 2 
are 0(L _1 ). Consequently, they may be neglected. The 
additive constant (the term proportional to c) can be 
neglected as well. With this in mind we determine: 

UH kin U ] =<52) 
J2{--(P u P ) 2 --+--(Pp) 2 --}+Wn + 2 PlPn ■ 



+ / dxdx 1 g(x — x')p\ J {x)pYi{x') + irvph 



TTV F / dx 



The expression in the curly brackets has to be trans- 
formed: 



dx {:{plf: +:(plf: } = j dx {: p%: +2: (pj) 2 : } . (53) 
The second term can be written as: 

2 J dx:(pl) 2 := (54) 

2 j dxdx'dx"v(x — x')v(x — x") : p p (x')p p (x"): = 
2 J dx'dx"{v*v)(x' -x") :pp(x')p p (x"): = 
2cl J dx:p 2 p :+h£l 



(4) 
int 

= H TL + nv F h{nl, 
v F = %(1 + 2c v ) = %(1 + 2v 2 ). (62) 

It is clear from eq.(61) that up to the linear combina- 
tion of irrelevant operators h-^. the classical Tomonaga- 
Luttinger model is unitary equivalent to the model of free 
fermions with the Hamiltonian H k i n - This was shown in 
[9]. 

At first sight it appears that our ability to map the 
Tomonaga-Luttinger Hamiltonian on the free fermionic 
quasiparticles is incompatible with the known fact that 
there is no quasiparticle pole in the Tomonaga-Luttinger 
single-electron propagator. The resolution of this para- 
dox is simple [9]: the overlap of the physical electron 
state and the quasiparticle state is zero. Thus, the singlc- 
quasiparticle propagator 

G qp (x, t) = - (0| % t)4(0, 0)} |0) , (63) 

V> p (x,t) = c tH ^iP p {x)c- tH ^, (64) 

does have a pole but such an object has no physical sig- 
nificance. Because of the orthogonality catastrophe the 
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physically important single-electron propagator 

Gtl(x,t) = -(0 TL |r T {^(x,T)^(o,0)}|0 TL ), (65) 



-T 



, (66) 



|0 T l) = U\0)- H tl ground state, (67) 

is not connected to the quasiparticle propagator in a 
Fermi liquid manner: £tl(<7,^) ZQ qp (q,u;)+Q reg (q,Lu). 
The Tomonaga-Luttingcr single-electron propagator was 
derived with the help of representation (61) in Ref.[9]. 

It is interesting to establish the relation between the 
Tomonaga-Luttinger Hamiltonian parameters vf, g and 
the free fermion parameter %. If we define go = 9q\q=a = 
4ttvfu v we can write: 



vf 




( .90 



\2lTV$ 



(68) 
(69) 



The latter formula is very well known in the bosonization 
theory. It gives the dependence of the plasmon velocity 
on the bare Fermi velocity vf and the interaction con- 
stant go- The latter interpretation does not contradict to 
eq. (69) for the plasmon mode in the fermion system with 
the Hamiltonian if kin propagates with the velocity vf- 

Our next task is to find UH n \W. We prove in Ap- 
pendix that: 



:(V^)(V^):-iv 2 p P 



4tt 2 



(70) 



This expression may be familiar to those who work with 
the bosonization: it is the bosonic form of fi nl [1, 11]. 
As a consequence of eq.(70) one derives: 



UH al U^ = ^ £ / dxU: Pp- Uf - 

The expression under the integration sign is transformed 
as follows: 

U:p 3 p {x):U*= (72) 
lim U {: p 2 Jx): p p (x') - 2b p {x - x')p p (x)\ C/ f . 

Using our result for U: p 2 :W we derive: 

U{:p 2 p :p' p -2b p (x-x')p p }U^ =(73) 

[: {p u P + p% + 5p p0 } 2 -- +ca 2 ] {( P ;y + { P \y + s Pp0 } 

-2b p (x - x') {pi + p"_ p + 5 Pp0 } . 

In this equality p' p = p p (x'), = p%> v (x'). 

The cubic term in eq.(73) has to be normal ordered: 

■■{P U P + P% + 5p P o} 2 :{(p u p y + (P V - P Y + S Pp0 } = (74) 

: {p u p + p\ + s Pp0 } 2 {(p u P y + ( P %y + 6p p0 y. 

+2[(u * b p * u)(x — x) + (v * 6_ p * v){x — x')] x 

{Pp+P- P + Sp P o}- 



Next one has to substitute eq.(74) into eq.(73) and collect 
similar terms. Using eq.(47) and definition eq.(51) it is 
easy to show now that: 

U:p 3 p .U^ = (75) 
: {P U P + P-p + S Pp0 } 3 : +3cA 2 {p u p + p% + 5p p0 }. 

Therefore, up to terms 0(L°) we find: 



£ J dx(:{p; + p%} 3 :+3cA 2 p p )(76) 



Now we need to expand the third order binomial: 

E : K + p-p} 3 -- = E [= (p u p) 3 -- +■ (p v P ) 3 -- ( 7? ) 

p p 

+3:(p;) 2 :p% + 3:( P ;) 2 :p u _ p ], 

and process the terms one by one. The first term: 

:{plf:=:{p p + p w p f:= (78) 
^:+3:p^:+3:p p ( P p 2 :+:(p^ 3 :. 

We use the same trick which was exploited while obtain- 
ing eq.(54): 

J ' dx:{p u p f: = (1 + wof J dx:p 3 p : + (79) 

J dx[{3:p 2 p p p ":-3w :p 3 p :} + 

|3:p p (^) 2 : -3w 2 : p 3 p . } + {: (p w p f: -w 3 :p 3 p . }] . 

As with eq.(54) it is possible to convince ourselves that 
the expressions in curly brackets have higher scaling di- 
mension than the scaling dimension of :p 3 :. Indeed: 



■P 2 p Pp--^' 



:p 3 p .= J2<^ 2n --P 2 p ^ 2n P P --- 



(80) 



71=1 



A 



In 



(2n)! 



J dxx 2n w(x). (81) 



Dimension counting shows that the operators under the 
summation sign have the scaling dimensions of 5 and 
higher. This is bigger than 3, the dimension of :p p :. 

Acting in a similar fashion we establish for the second, 
third and fourth terms of eq.(77) that: 



J dx:{pl) 3 : = v 3 j dx:p 3 p . + J dx[.{ Pp \) 3 : -v 3 :p 3 p .} ,(82) 
J dx(:(p;) 2 :pl p +:(p;) 2 :pl p )=(83) 
(u 2 v +v 2 u ) J dx: p 2 p . p_ p + 

j dx ( ■{Pp) 2 -P- p - u l v o : P 2 p -P- P + 
: (filf-P-p - v l u o-P 2 P -P- P ) , 
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where the expressions in square brackets has the scaling 
dimension 5 or higher. 

Combining the above results we have: 



+ 47T 2 t4 „ x 



dx(: P l+:pi:)+ (84) 



27T ffex / dx [: p L : p R +: p R : p L j + —^—h\ n l , 



«f,c X -("o + «o)«f, (85) 
<4 = 2 ™f(^o+^o)- (86) 

The superscript in the notation h^ t reminds us of the 

(5) 

fact that h\ n l is a series of operators whose lowest scaling 
dimension is 5. 

Further, combining eq.(43) and eq.(70) we write: 



/A 2 ~/ 
dx V$W p : +i<4P (: ^W p ): -: (V^)^ p : ) P- P } + (hg> + 3cA 2 (iV L + iV R )) .(87) 

p 

Therefore, using the above formula and eq.(61) one finds the following expression for H ex = UH U^: 

H cx = ivp J dx (: ^ h S7^ h : -: V>£VV> R : ) + J dxdx 1 g(x - x')p L (x)p R (x') + (88) 

<,ex^ / dx:ViPlViP p :+ig' ex J2P [ dx ( : V> P (W> P ): -: {W P )%- ) P- P + M^l + N K ) + mfchf 
P J P J 



'' ! _i_ ^F fa (5) 



Since at v F ^ the particle-hole symmetry is broken the 
bare value of the chemical potential is shifted: in the for- 
mula above p = 4tt 2 cv f A 2 whereas p = in Hamiltonian 
(8). 

It is important to realize that v' F ex , g' ex and g in eq.(88) 
are not independent. By dividing eq.(85) on eq.(86) we 
obtain: 



J F,cx 



2tt 



l 



"0 



Since 



we establish: 



u 2 + Vq 

Mo «0 



U V 



A-KVv, 



v F,cx — 



2v F 

go 



90 



1 

2^ 



(89) 



(90) 



(91) 



The deduced results can be summarized in the from of a 
theorem. 

Theorem 1. Consider the Tomonaga-Luttinger model 
with nonlinear dispersion and the irrelevant operators 
h-*f^ whose scaling dimensions are 4 and higher. Such 
model is unitary equivalent to a system of free fermions 
with Hamiltonian (8) provided that parameters v' F ex , g' ex , 
v F and g satisfy eq.(91). 

The theorem trivially implies that the Tomonaga- 
Luttinger model with Hamiltonian (88) whose parame- 
ters satisfy relation (91) is exactly soluble. This allows 
one to find the density-density propagator of eq.(88). 
Since the ground state |0tl) of Hamiltonian (88) is given 



by eq.(67) we derive: 

(0tl| Wt) + pr 9 (t)}{ Pl _ 9 (0) + p R - q (0)} |0tl) K92) 
(0| UH PLq (T) + PRq (T)}UUHpL- q (0) + P R- q (0)}U |0) . 

Here, by the Green's function definition, p pg (r) = 
exp(ri? ex )p P9 exp(— tH cx ), that is the evolution is con- 
trolled by Hex, eq.(88). The action of W = If- 1 on the 
density and evolution operators is: 

U^(e rH " Ppq e- TH ")U = e TH °(u qPpq - v q p_ pq )e- TH °(93) 

Note the minus sign in the above transformation law. It 
differs from eq.(23) since here we apply W rather than 
U. Observe as well that the evolution of the density 
operators on the right-hand side of the above equation is 
set by the Hamiltonian (8). Thus: 

(0tl| {pl 9 (t) + PRq (r)}{p h - q (0) + p R - q (0)} |0 TL ) =<94) 
(«, - v q f (0| { PLq + PKq ]e- THa {p^ q + p K _ q } |0) . 

The latter expectation value is easy to find. It is just a 
product of two free single-quasiparticle propagators cor- 
responding to Hamiltonian (8). Therefore, Matsubara 
propagator for the Tomonaga-Luttinger model (88) is 
equal to: 



V^(r) = (u q -v q ) 2 V° q (r), 



(95) 



where T>° is the propagator for the Hamiltonian H spec- 
ified by eq.(13). At low |g| this formula may be cast in a 
traditional for the bosonization literature way: 



K= (u - w ) 2 = 



2ttv f - ff 
27TWF + go ' 



(96) 
(97) 
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Here /C is the usual Tomonaga-Luttinger liquid parame- 
ter. 

The retarded propagator D q * and the spectral density 
B£ = -2ImD- are: 

p CX g f (£ Fg - £ Fg 2 ) 2 - + *Q) 2 ] f9fn 

^ 47ri^g S 1 (w F g + « F <7 2 ) 2 - (w + iO) 2 J ' k ' 

B ?« = ^- q {0(" 2 -^-^ 2 )- (") 

$ (w 2 - (5 F <7 + «f<? 2 ) 2 ) } s S n w - 

Due to exact solubility the spectral density is zero ev- 
erywhere outside the interval v F \q\ — \v' F \q 2 < \u)\ < 
i F \q\ + \v' F \q 2 . 



III. NON-TECHNICAL OVERVIEW OF THE 
DERIVATIONS 

The presentation of the previous Section was rather 
technical. Below we will give a less formal account of the 
derivations performed so far. 

It is obvious that the major source of technical compli- 
cations are the zero modes and the irrelevant operators 
/j(4,5) _ -j- ne same time, they are the least important 
parts of the derived Hamiltonian (88). Thus, it is useful 
to redo the calculations of the Section II paying no atten- 
tion to the zero modes and highly irrelevant operators. 
That way we will be able to capture easily the most im- 
portant features of the algebraic manipulations without 
losing too much significant information. 

Since the material presented in this Section is a refash- 
ioning of the discussion given above the reader should be 
prepared for some degree of redundancy. 

Our technical goal is to start with the Hamiltonian of 
the free fermions with nonlinear dispersion Hq, eq.(8), 
and transform Hq into H ex , eq.(88). We begin the exe- 
cution of this program by writing H n as: 



f i - 2/ \ 4:TT 2 V V ,. , 

H = dx^nv F :p p {x):+^—:p p {x):. 



(100) 



The validity of this equation can be checked with the 
help of eq.(43) and eq.(70). We already pointed out in the 
previous Section that the above equation may be thought 
of as the bosonic form of H n [1, 11]. 

As in Sect. II we want to transform Ho with the help 
of U and find UHqW . Operator U acts as a Bogoliubov 
rotation on the density operators, thus we may write (see 
eq.(35), eq.(50) and eq.(75)): 

Up p (x)U f = u Q p p (x) + v p- p (x) + (101) 
U:p 2 p {x): C/t =:(u p p (x) + v p- p {x)) 2 : + . . . , (102) 
U:p 3 p (x): C/ f =:(u Q pp(x) + v Q p- p (x)f: + . . . , (103) 

where the ellipses stand for (a) additive constants, (b) 
zero modes terms (c) irrelevant operators with high scal- 
ing dimensions and (d) operators which are reduced to 



zero modes upon integration over x. As it was explained 
above such terms do nothing but cloud our view. Thus, 
we will neglect them. 

Applying the rules eq.(101)-(103) to eq.(100) the fol- 
lowing is derived: 

UH U^ = dxJ2^F'-( u oPp + voP- P ) 2 - (104) 
J p 

47T 2 W F ., 

H g — :(uoPp + vop-p) : + ... 

= [ dxJ2^ F (u 2 + v 2 ):p 2 p .+^-^(u 3 + v%):p p : 
J p 

+2irv F u v PpP-p + 4tt 2 v f (uIv + v 2 u ):p 2 : p_ p + 

where the ellipsis stand for the zero modes and highly 
irrelevant operators. One notes that the term :p 2 : is pro- 
portional to ip-.ip^ Vip : and, thus, the coefficient in front 
of :p 2 : becomes v F : v F — vf(uI + v 2 ). Likewise, the co- 
efficient of pl Pr becomes go (see eq.(60)), the coefficient 
in front of :p p : p_ p is <? ex (see eq.(86)) and the coefficient 
of :p p : is v' F ex (see eq.(85)). Consequently, we obtain: 

H cx = UH U^ =(105) 

dx ipv F :4>lVip p : +v' F cy /^7ip p Vip p . +^9oP P P- P 
p 

+*P3L (■ ^(VV P ): -: (W> P )Vy ) P- P + ■ ■ ■ ■ 

The previous equation coincide with eq.(88) up to zero 
modes and highly irrelevant operators. The coefficients 
g' c ^ and v' F satisfy eq.(91). This concludes our nontech- 
nical overview of the exactly soluble Hamiltonian deriva- 
tion. Transcending many technical details discussed in 
Section II we discover that the method of calculating i? ox 
is in fact fairly simple. 

Another question we want to discuss here is the mean- 
ing of the exact solubility line eq.(91). Consider the 
generic Hamiltonian Hq: 



I 



Ho = Htl + H n \ + H in 



(106) 



whose parameters do not satisfy any specific relation. 
It is possible to find the transformation V of the form 
eq.(14) which diagonalizes the first term of the above 
equation. For the detailed calculations one should con- 
sult Rcf.[9] or our derivation of eq.(61). The heuristic ex- 
planation, however, is easy. Consider the 'bosonic' form 
of H TL : 

Htl^^YI [PpiPp-1 + ^~ p P" p -P-i) + ' • ( 10? ) 

pq ^ F ' 

Application of V produces the following: 

VHtlV* 1 = (108) 

7TV F 

T 



7TV F X 



pq 



u 2 q + v 2 + ^ 



^L-(u 2 q + v 2 ) + 2u q v q 



PpqPp-q + 

PpqP—p—q + • • • 
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By adjusting u q and v q we can kill the interaction term 
PpP-p- For this to happens u, v must satisfy the equation: 



2irv F 
or, equivalently: 



(u q +V 2 q ) + 2UgV q = 0, 



th 2a„ 



9q 



(109) 



(110) 



q 2irv F ' 
If this condition is met we have: 

VH TL V^ = iv F J dx (: V>£VV L : -: <V^ R : ) + . . (711) 

where the product p p p p was expressed in terms of the 
fermionic operators ip,^ with the help of eq.(43). 

There exists another transformation V' of the general 
form (14) but with a different set of coefficients a' which 
diagonalizes the sum (H n i + H[ nt ): 

V'(H nl + H' int )(V'y =v' F f <ix]r:V^VVy + ..(112) 

J v 

For this to take place the parameters a' must satisfy the 
relation: 



2-Kv'p + g'(l + 2cth 2a' ) = 0, 



(113) 



where a' equals to a' q at q = 0. The derivations of the 
two above equations is similar to the proof of equations 
(110) and (111). The idea of these derivations is to write 
V'(H nl + H' int )(V'Y in terms of p's: 



V'(H nl + H! nt )(Vy = V 



I 



(:p£:+:^:)(114) 



+2ng'(:plp R +:p 2 R :p L ) + ...\(V r ) 

and adjust a' Q of V' in such a way as to remove : p 2 p : p_ p 
terms. This requirement lead us to eq.(113). Once the 
latter terms are removed what remains is proportional to 
: : + : Pr '■■ This expression is proportional to H n i 
and, therefore, is diagonal. 

If the values of v F , v' F , g and g' are arbitrary then 
V ^ V . Consequently, the whole Hamiltonian Hq can- 
not be diagonalizcd with the help of the canonical trans- 
formation of the form eq.(14). However, if we demand 
that V = V or, equivalently, a = a' such diagonaliza- 
tion becomes possible and the system becomes soluble. 
One can check that the exact solubility condition (91) 
may be derived from the equation a = a' . Thus, one 
can say that the exact solubility condition guarantees the 
simultaneous diagonalization of the whole Hamiltonian. 



IV. GENERIC TOMONAGA-LUTTINGER 
MODEL WITH NONLINEAR DISPERSION 

In this Section we will discuss how the above findings 
can be applied to a generic model of the one-dimensional 
spinless fermions. 



The most general Hamiltonian is: 

H G = H Th + H nl + H! mt + H (4) , (115) 

where stands for operators whose scaling dimensions 
are 4 and higher. Here no special relation between v' F and 
g' is assumed. Neither do we suppose that H^> is equal 
to the linear combination of and M 5 -* from eq. (88). 
To avoid clutter we do not show the chemical potential 
term explicitly for it can be accounted trivially. 

It is convenient to write Hamiltonian Hq in the form: 



H G = H ex + SH nl + SH^, 
6H nl = Sv' F J dx ^:VV>tVVy , 



Sv' F 



(116) 
(117) 

(118) 



The term 5H^ contains those operators whose scaling 
dimension are 4 and higher and which are not absorbed 
into H cx . The constant 8v' F and SH^ measure the de- 
viation of the Hamiltonian Hq from the exact solubil- 
ity line. We want to apply the transformation V of the 
form eq.(14) to Hq such that the transformed Hamilto- 
nian contains no marginal interaction operator H mt . The 
remaining irrelevant operators can be treated perturba- 
tively. 

We already know how the first two terms of eq.(116) 
are transformed under the action of V. The unknown 
entity is VSH^Vl To discuss the action of V on 5H^ 
it is useful to introduce the notation: 

fckin =ijdx{-. VJVVl: 4W> R : ) , (119) 

hint = J dxp L p R , (120) 

h[ nt = J ^^w_ p (:^ (V ^):-:(V^)V p :),(121) 

'dx]T:V^V^:.(122) 



These operators exhaust the list of operators whose scal- 
ing dimension is lower than 4 and we can write the most 
general expression: 

VSHWV* = (123) 

Adn^kin + Ant^int + Pn\h n l + A'nt^int + SH^ . 

The coefficients /3's are functions of a's which specify V. 

Now we must consider two cases: (i) when (5- lrit = 
and (ii) when (3i nt ^ 0. 

Let us start with (i). The transformation V with 



a 



9o 



-Ith- , 
2 V 2irv F 



(124) 



acts on Hq as follows: 

VHqV^ =H + (Sv' F h rd + 5~g'h[ nt ) + VSH^V^ (125) 
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= [H + Adnft-kin + {5v' F + /3„l)/l„l] 

+(8g' + p[ nt )hi nt +Sm 4 \ 

Sv'p = (u 3 - vl)6v' F , (126) 
Sg' = 2tt(vqU - u 2 ,v )5v F . (127) 

In the first line of eq.(125) the expression in the round 
brackets corresponds to V8v' F h n \V^ . 

Eq.(125) can be reformulated in a more compact way: 

VH G V ] = Ho + 5H irr , (128) 
H Q =H Q + /?kin/ikin + (H + Pm)h nh (129) 
5H lir = (5~g> + 0{ nt )ti int + SH™. (130) 

Hamiltonian Ho is quadratic in the fcrmionic field opera- 
tors. It is similar to Hq but has renormalizcd parameter 
values. Likewise, SH^ is the renormalized version of 
8H^. 

We see that deviation from the exact solubility line 
eq.(91) in Sv' F has two effects: (a) renormalization of v F 
and v' F and (b) introduction of the interactions between 
the quasiparticles (operator <57J; rr ). Fortunately these 
interactions are irrelevant in the renormalization group 
sense. Thus, their influence could be accounted pertur- 
batively. For that reason let us introduce the notation: 

H qp = VH G V\ (131) 

which suggests that the transformed operator VHqV^ 
represents the Hamiltonian of the weakly interacting 
quasiparticles. 

Next we study the case (it). This situation is not much 
more difficult than (i). From the field-theoretical point of 
view /3j n t ^ means that irrelevant operators contribute 
to the renormalization of the marginal operator coupling 
constant g. The correct way to address the problem is 
to fix ao by a generalized version of eq.(124). The value 
of ao must be determined by the requirement that the 
quasiparticle Hamiltonian 

H qp = 7kin^kin + 7intftint + 7nl^-nl + 7int^int ( 132 ) 

+5H^ 

(7's are functions of ao) contains no /i; nt : 

7intM=0. (133) 

It is easy to demonstrate that eq.(133) is equivalent to 
eq.(124) if Ant = 0. 

Once ao is determined according to eq.(133) the ex- 
pression (132) attains the form of eq.(128). Thus, we 
prove 

Theorem 2. A generic Tomonaga-Luttinger model 
with nonlinear dispersion is unitary equivalent to a sys- 
tem of the free quasiparticles with irrelevant interactions. 

V. CALCULATIONS OF OBSERVABLES 

What does Theorem 2 imply for evaluation of the phys- 
ical properties of HqI The irrelevance of corrections 



8H lvv means that the ground state properties of Hamilto- 
nian Hq can be approximated by the ground state prop- 
erties of Ho- Should we be interested in corrections those 
could be found perturbatively. 

Further, for an observable 0(x) the correlation func- 
tion 

Qo = - (0 G | % {0(x, T)0t(o, 0)} |0 G > , (134) 
0(x,t) =c THc -0(x)e~ THG , (135) 
|0 G ) - H G ground state, (136) 

may be expressed as: 

Go = - (0qp| % {o(x,T)Ot(0,0)}|0 qp ), (137) 

6{x,t) = e rH ™ (V0(a;)V + ) c- tH ^, (138) 
|0 qp ) — ground state of H qp . (139) 

If the observable O is such that VOV^ has a very com- 
plicated form the practical evaluation of eq.(137) may be 
difficult to carry out. An example of such situation is 
O = tp p . At present the author does not know how the 
single-fermion propagator for the case of nonzero v' F can 
be calculated. 

On the other hand, if the object VOV' has some sim- 
ple representation in terms of the quasiparticle degrees of 
freedom one can evaluate Qo with the help of the usual 
perturbative expansion in orders of 5H m . For example, 
consider calculation of the density-density propagator: 
O = (pl + Pr)- In that instance VOV^ is a linear com- 
bination of the density operators. The zeroth order ex- 
pression for such propagator is: 

V qu} =V q l + 0((5H m f). (140) 

The corrections in orders of 5H{„ could be determined 
with the help of field-theoretical techniques. 

Note, that the actual calculation of 0((8H- m ) 2 ) cor- 
rections may be quite nontrivial. Indeed, numerical ev- 
idence [14] and heuristic analytical arguments [15] indi- 
cates that D quJ of the generic model diverges stronger 
than D q *. This suggests that the perturbation theory 
for D qLU requires a reliable resummation technique. 

However, we are not always interested in the prop- 
agator itself. In fact, there are situations when a few 
low-order diagrams contributing to D quJ would suffice. 
For example, consider some quantity r proportional 
to an integral over D qu} . [An important example of 
such quantity is the Coulomb drag resistivity [6]: r cx 
J dqdwf(q, w)(Im D quj ) 2 , where / is some temperature- 
dependent kernel.] Unlike the propagator itself, r re- 
mains finite. This is so because the divergences of D qu 
due to 0((SHi rl ) 2 ) corrections are very weak. Thus, the 
integration of D quJ can be successfully performed with 
the result: 

r = r cx + 5r, (141) 

where quasiparticle interaction correction \Sr\ < 00 may 
be non-analytic. 
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Regarding the quasiparticle interaction corrections due 
to (>-ffi rr we would like to make a remark. Operator 
(5g' + (3{ nt )h' int is the least irrelevant part of 8H m . Con- 
sequently, we expect that at small momenta the correc- 
tions due to (Sg' + f3' int )h[ Dt are stronger than those due 
to SH^ 4 \ Therefore, the most important contribution to 
the low momenta corrections are cx (Sg' + (3[ nt ) 2 . Thus, 
there is a hypersurface in the generic one-dimensional 
model parameter space where the low momenta correc- 
tions drastically weakens. This hypersurface is fixed by 
the equation 

Sg + A'„t - 0. (142) 

We can interpret this as a manifestation of the exact 
solubility line eq.(91). 

As a practical application one can calculate the hyper- 
surface (142) for the spinless Hubbard-like model 

■ffhub = ^2t(i- j)c\cj +^^n 4 + u^n 4 n l+ i(143) 

i,j i i 

in the limit of small interaction u<t. Using the decom- 
position: 

Ci = Mx t )e- lkFX > + ip R { Xi )e ikFXi , (144) 

one can calculate vf, v' f , g and g' in terms of t(i) and 
u. As soon as this task is done the equations for those 
quantities must be substituted into eq.(91) and the limit 
u — > is taken. This gives: 

v F,ex = a^Fctg (fcpa), (145) 

where a is the lattice constant and the values of wp and 
v' F ex are functionals of t(i). 

What are the phcnomcnological manifestations of 
eq.(145)? As explained above, on the surface fixed by 
the latter equation and in the limit it <C i the scattering 
of the quasiparticles is drastically reduced. As a con- 
sequence, the propagator D qix] calculated for the Hamil- 
tonian satisfying cq.(145) differs qualitatively from the 
propagator of the generic Hamiltonian eq.(143). For ex- 
ample, the divergence of the propagator due to quasipar- 
ticle scattering might be weaker. Such prediction may be 
tested numerically. 

VI. COMPARISON WITH OTHER WORK 

The issue of nonzero v' F has a long history. Tradi- 
tionally, it has been approached within the framework 
of the bosonization [1, 11, 13]. The disadvantage of the 
bosonization is that the perturbative expansion in orders 
of v' F is highly singular. 

To circumvent this latter difficulty the authors of 
Ref.[12] used the random phase approximation to find 
D qu} . Their spectral function differs qualitatively from 
eq.(99). We believe that such a result is an artefact of 
the random phase approximation. 



Alternatively, in several recent papers [6, 7] the use of 
the 'fermionic' methods was advocated. The idea behind 
them is that away from the mass surface the perturbation 
theory in g is permissible without bosonization. That is, 
the authors chose to work with the original fcrmions but 
at a price: the perturbation theory results cannot be 
applied straightforwardly near the mass surface. 

The proposed method was used to calculate the high- 
frequency tail of the propagator. The latter technique 
is complimentary to the exact solution: it is easily im- 
plemented at high frequency but it is problematic at the 
mass surface where the exact solution works. 

In Ref. [15] the heuristic argument was put forward 
suggesting that the spectral function B of Hq diverges 
at some line of (q, u>) plane. This finding is corroborated 
qualitatively by the numerical work [14]. How does this 
compare against calculations presented here? 

Obviously, our spectral function B cx of the exactly 
soluble model is finite. This docs not contradict to 
Ref. [14, 15]: the propagator of the generic model Hq is 
affected by the quasiparticle interaction which induces 
the extra divergence. Thus, the propagator calcula- 
tion requires an adequate resummation technique, as we 
pointed out in the previous Section. We may hypothe- 
size that the argumentation of Ref. [15] may be directly 
applied to the quasiparticle Hamiltonian H qp . Such 'mar- 
riage' of the two approaches possesses a clear advantage: 
both marginal and irrelevant interactions are accounted 
for to all orders. Further speculations are to be post- 
poned until thorough investigation. 

To conclude, we have shown that the Tomonaga- 
Luttinger model with specific parameter values is equiv- 
alent to the free fermions. At this point the Tomonaga- 
Luttinger Hamiltonian is exactly soluble. Away from the 
exact solubility the model may be mapped on a system of 
the fermions with irrelevant interactions. Although, un- 
der such circumstances the model is not soluble exactly, 
yet, it can be investigated with the help of the perturba- 
tion theory. 
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VIII. APPENDIX 

In this Appendix we prove certain relations required 
by the presentation in the main text. 

First, we derive eq.(39). Our initial step is to prove: 

■Pp(x)p P (y)- = ■ipt(x)i) p (x)ipl(y)i) p (y): + 
s p (x - y):ip p (x)ipl(y): +s p (x - y):ipl(x)ip p (y): . (146) 
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The proof goes as follows: 

Pp{x)p P {y) = ■ ^ p {x)%(x): :^ p (y)%{y): = (147) 

4( x ~ y)+-i ) li x )%{ x )i ) l{y)%{y)- + 

s p (x - y):tp p (x)tpl(y): +s p (x - y):ip P (x)ip p (y): , 

where we used eq. (10). Taking into consideration 
eq.(41) and b p = s p one obtains the desired formula 
eq.(146). 

Once eq.(146) is established eq.(39) follows. We send 
y — > x and note that :ip P (x)'tp p (x)'tpj ) (y)^p p (y): = if x = 
y. Indeed: 

■4l{x)%{x)^ p {y)%{y):= (148) 
-^t(y)^ P (x)^l(x)^ p (y): ■ 

Since the quartic term vanishes we need to evaluate the 



quadratic terms of eq.(146) only. To calculate them at 
x = y one can utilize the L'Hospital's rule: 

lim [ Sp (x-y){:i> p (x)i>l(y):+:iPl(xW p (y):}] =(149) 



ip 
2^ 



{:i> p (x)^l(x):+:i>l(x)ViP p (x): } , 



where we took into account that 



d 



dy \s p (x - y) 



2tt 

ip 



(150) 



It is obvious now that eq.(39) is valid. 

Next, we want to prove eq.(70). The procedure is sim- 
ilar to the one we have used just now. It is possible to 
show: 



J 



: Pp{x)Pp{y)p P {z): =: %{x)ip p (x)i>l(y)ip p (y)ipl(z)ip p (z): (151) 
+s P (x - y) {: ipl(x)ip p (y)ipl(z)ip p (z): -: ipl(y)ip p (x)ipl(z)ip p (z): } 
+s p (x - z) {:ipl(x)ip p (z)ipl(y)ip p (y): -: ipl(z)ip p (x)ipl(y)ip p (y): } 
+s p (y - z) {•i>l(y)i) p {z)i>lix)'il) p {x)\ -: ipl{z)il) p {y)^l{x)i} p {x): } 
+s p (y - z)s p (x - z){:i) p (x)i>l(y): -:il>l(x)ip p (y): } + s p (x - y)s p (x - z) {:ip p (y)ipl(z): -:il>l(y)i> p {z): } 

-s p {x - y)s p (y - z) {:il) p {x)i>l{z): -: ipl(x)i> p (z): } . 

I 



Like eq.(146), this equality was deduced by writing ppp as 
a product of six fermionic fields which is then expressed 
as a sum of the normal ordered expressions. If we put in 
the above equation x = y = z only terms quadratic in 
the fermionic fields would survive. This is a consequence 
of the fermionic field anticommutativity (see eq.(148)). 
Since we are interested in x — y — z limit we have to 
handle the quadratic part of eq.(151) only. 
Let us define an operator-valued function: 

G(x,y,z)= (152) 

s P (y - z)s p (x - z) {:ip p (x)il>l(y): -:^l{x)^ p {y): } 

+s p (x - y)s p (x - z) {: iP p (y)iPl{z): -: i> P {y)^ p {z)-. } 

-s p (x - y)s p (y - z) {:ip p (x)ip P (z): -:ipl(x)ip p (z): } . 

It coincides with the quadratic part of eq.(151). There- 
fore, : Pp(x): = G(x, x, x). We will demonstrate that 

G(x, x, x) = :pp(x): = (153) 

A ^ (V^ (x) )(v^(x)):-iv 2 p P (x)) . 

We begin the proof by an observation: 

s p (x - y)s p (y - z) = s p (x - y)s p (x - z) + (154) 
s P (y - z)s p (x - z). 



The formula can be checked directly with the help of 
eq.(ll). We transform eq.(152): 

G(x,y,z)= (155) 
s p (x - y)s p (x - z) {: [^(x) - V>p(j/)] %{z): +h.c.} - 
s p (y - z)s p (x -z){: [^(y) - 4>l(z)] ip p {x): +h.c.} . 

It is convenient to define another operator-valued func- 
tion F: 

G(x, y, z) = s p {x - z) [F(x, y, z) - F(y, z, x)] , (156) 

F(x,y,z)= (157) 
s P (x -y){- [i>l(x) - i>l(yj\ ip p {z): +h.c.} . 

Function is well-defined for any values of x, y and z. 

It is important to note that F is invariant under the 
transposition of its two first arguments: 

F(x, y, z) = (6{x -y)- s p {y - x)) x (158) 
{-: [^(y) - ^(z)] ip p ( z ); +h.c.} = F(y,x,z), 

since 8(x — y):[ip P {x) — V|(y)]V'p( z ) : = 0. Consequently, 
we can rewrite the expression for G: 

G(x, y, z) = s p (x - z) [F(x, y, z) - F(z, y, x)] . (159) 
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The expression in the square brackets vanishes at x = z. 
Thus, function G is well-defined everywhere. Sending 
z — > x we obtain with the help of the L'Hospital's rule: 

G(x,x,x) = ^— [F(x, x, z) — F(z, x, x)] \ z=x . (160) 

Applying the L'Hospital's rule the following equality is 
derived: 

F(x,x,z) = £-.{: {V1>l(x))%(z): +h.c} . (161) 
Consequently: 

§-F(x,x, z)\ z=x = £: (V^( x) ) ( v^(x)): . (162) 



We calculate: 

^F(z,x,x)\ x=x = £- |:(V 2 ^ (a; ))^( x) ; +h . c .}(i 6 3) 
= £~ {^PM - 2:(V^ (x) )(v^(x)): } , 

by applying the L'Hospital's rule yet again. 

Combining the results for G and derivatives of F we 
prove eq.(153) which, in turn, proves eq.(70). 
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